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Abstract 

This paper describes an approach based on using both 
bias momentum and multiple control moment gyros for 
controlling the attitude of statically unstable thrust- 
levitated vehicles in hover or slow translation. The sta- 
bilization approach described in this paper uses these 
internal angular momentum transfer devices for stabil- 
ity, augmented by thrust vectoring for trim and other 
“outer loop” control functions, including CMG stabi- 
lization/desaturation under persistent externanl distur- 
bances. Simulation results show the feasibility of (1) 
improved vehicle performance beyond bias momentum 
assisted vector thrusting control, and (2) using control 
moment gyros to significantly reduce the external torque 
required from the vector thrusting machinery. 

1 Introduction 

In this paper, we focus on the vehicle controls aspect 
of the Vertical Take-off and Landing (VTOL) capable 
NASA Flying Test Platform (NFTP) as reported earlier 
in [1], This earlier paper describes the recent interest 
in flight vehicles, capable of VTOL and hover, that are 
statically unstable or neutrally stable. In particular, it 
discusses two specific difficulties inherent in the state of 
the art approach to vehicle control by thrust vectoring 
and differential thrusting, possibly with additional aero- 
dynamic features implemented to enhance stability in the 
presence of gusts. For convenience, we reiterate these con- 
cerns. First, because of the airframe’s lack of stability, the 
controller must operate with high authority and band- 
width in order to contain and suppress excursions in at- 
titude. Accurate system models are generally required to 
successfully design such controllers. The second difficulty 
is that such models are not available. This is because 
the active thrust vector variations required for stabiliza- 
tion in the presence of disturbances induce unsteady flow 
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phenomena that distort the control command in a man- 
ner that is currently prohibitively difficult to model for 
control design. These difficulties come together when the 
control engineer attempts to design the required high- 
performance controller while ignoring or under-modelling 
the effect of the unsteady flow phenomena on the closed- 
loop dynamics. These difficulties become more evident 
for VTOL vehicles with smaller inertias because they re- 
spond more rapidly to disturbances so that correspond- 
ingly more rapid feedback control is necessary for vehicle 
stabilization. 

The earlier paper [1] describes a novel approach to 
circumventing the problem of control-induced unsteady 
aerodynamics on the closed loop stability during vehicle 
stabilization. A tradeoff is described in which a significant 
level of stored bias momentum in the 2 -body axis pro- 
vides directional stability and disturbance response mit- 
igation at the cost of additional wheel weight and com- 
plexity of the dual spin dynamics for the system. The 
sizing of the bias momentum level for given nominal ve- 
hicle inertias, disturbance levels, and accuracy levels have 
been described. Detailed nonlinear simulations based on 
closed loop controlled response under strong turbulent 
winds and significant in-flight payload variations clearly 
demonstrated the significant advantages in terms of atti- 
tude and trajectory robustness. 

In this paper, we demonstrate the feasibility of using 
a set of variable speed control moment gyros (VSCMG) 
to obtain controlled performance beyond what is possible 
through a bias momentum stability-enhanced system. We 
discuss a methodology for synthesis of a nonlinear con- 
troller for the system which incorporates vector thrusting 
and differential throttling, bias momentum wheel, and 
variable speed CMGs. Closed loop stability of the non- 
linear system is enforced by well known Lyapunov stabil- 
ity conditions. Control performances are described and 
compared to illustrate the improvements possible through 
the use of CMGs, especially in its ability to reduce the 
commanded torque to the vector thrusting and differen- 
tial thrusting machinery. The controlled performance of 
interest in this study is the attitude robustness of the ve- 
hicle during low airspeed or hovering operations during 
strong turbulent winds. 


1.1 Outline of Paper 

This paper complements the results reported recently in 
[1], In the next section, a research vehicle under devel- 
opment at NASA Langley Research Center for maturing 
this control technology is described. In addition, the dy- 
namical equations of motion about trim used specifically 
for control law design, and a flight dynamics and control 
framework are also given. Section 3 discusses controller 
synthesis for the system, which involves the integration 
of vector thrusting and CMGs for a dual spin stabilized 
VTOL air vehicle. A baseline set of nonlinear controllers, 
that can be readily designed, are described which satisfies 
the conditions for global asymptotic stability. In section 
4, we describe simulation results that compare the perfor- 
mance of a set of control laws that do not use VSCMGs 
to a second set that do. Section 5 provides conclusions. 

2 Control Design Model 

2.1 NFTP Vehicle Description and Mod- 
elling Assumptions 

The NFTP is designed as a laboratory-scale vehicle that is 
levitated by vectored thrust, additionally equipped with 
a suite of variable-speed CMGs and a momentum wheel 
that provides a substantial constant angular momentum 
vector aligned with the vehicle’s z-body axis. A CAD 
model is pictured as Figure 3, and dimensional and mass 
properties are given in [1], The salient features of the 
vehicle’s dynamics are: 

• mass roughly 14kg, maximum thrust roughly 200N 

• principal moments of inertia roughly (.7, .7, 1.4) Kg- 
m 2 

• z-body angular momentum bias roughly 18N-m-s 

• VSCMG set output torque capacity roughly lON-m 

• peak z-axis torque due to vanes are roughly 8N-m 

The net thrust vector is produced by the four ducted 
fans appearing in Figure 3, two of which rotate clock- 
wise, and two counterclockwise, in order to discourage 
yaw torque transients related to thrust variation, and to 
encourage a constant value for the overall z-body angu- 
lar momentum bias. Downstream from each fan is a pair 
of thrust- vectoring vanes. Net control force and torque 
vectors can be generated by differential thrusting and by 
deflecting the thrust vectoring vanes. These net control 
forces and torques are the result of adding trim control 
forces and torques and commanded thrust and torque 
deviations from trim, T v and t v . For simplicity, these 
thrust and torque deviations are directly commanded by 
the control system; a numerical procedure for efhently 
mapping these deviations to engine and vane settings has 


been developed, but is not reported here. In an actual 
application, these required thrust and torque deviations 
from trim must be implemented in such a way that will 
account for the effects of ducted fan and motor dynamic 
inertial coupling with the platform to which it is attached 
(see [3] for details). 

The momenta of the CMGs along their spin and gim- 
bal axes are, respectively, 


hw = t s 
he = i’ + Tg 


(1) 


where t s is the four-element vector of control motor 
torques on the wheels’ spin axes, and r g are the corre- 
sponding control motor torques on their gimbal axes. The 
vector tji contains angular momentum coupling from the 
platform and CMG wheel angular velocities. 

Finally, the z-body location of the center of mass is 
near, but above the plane formed by the points at which 
the four fans thrust vectors are applied. Because of this, 
the vehicle is slightly statically unstable. Reference [1] 
lists a number of potential applications of vehicles of this 
general type. A key feature of many of them is that the 
center of mass of the payload is above the fan plane, and 
in an uncertain and variable position in the x and y-body 
axes. 


2.2 Dynamics about trim 


Trimmed flight for the NFTP, in this paper, consists of 
maintaining zero airframe rotation rates and translational 
accelerations. The CMGs, in this case, have constant spin 
rates and gimbal angles; thus, the CMGs take no active 
part in trimming the vehicle. Trim is achieved by an 
appropriate selection of fan speeds, Cj s £ H ixl and vane 
angles, 9 £ lZ 4xl , an algorithm for which is described in 
[1]. Note that nonzero r„ is generally required in order 
to offset gravitational torque arising from the interplay 
between the attitude and the center of mass location. 

Denoting airframe translational and rotational veloc- 
ities as v and uo, and the CMG spin axis velocity and 
gimbal angle vectors as f 1 and ry, the dynamics of the 
system for excursions about trim can be written as 


MX = Z+U, 



( 2 ) 


where M and Z are complicated nonlinear functions of 
the rotational states and mass properties. A detailed 
presentation of (2) is under preparation [3]. The form 
of (2) is significant, because it permits the synthesis of 
the Lyapunov-based nonlinear controller described in the 
next section. 
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2.3 Vehicle attitude parameterization 


2.4 Control Framework 


With the target trim attitude chosen as a reference at- 
titude, the vehicle’s attitude denoted by body frame Fb, 
can be conveniently viewed as an excursion about a trim 
frame, F t , as shown in Figure 1. In particular, the in- 



Target (Trim) Frame 
Vehicle (Body) Frame 


i 

NED (Inertial) Frame 



Figure 1: Reference frames for attitude tracking. 

stantaneous attitude of the vehicle with respect to an 
inertial frame, F a , is viewed here as a successive rotation 
from inertial to trim frames followed by a rotation from 
trim to body frame. The direction cosine matrix for the 
former transformation is parameterized by Euler param- 
eters, (3 € TZ 4x1 , while the latter transformation is pa- 
rameterized using Modified Rodrigues Parameter (MRP) 
e € K 3xl [5], 


F t 

= C t0 (f3)F 0 

(3) 

F b 

= C bt (e)F t 

(4) 

F b 

= CF a ; C = Cbt(e)C to 0) 

(5) 


While the choice of parameterization of the directon 
cosine matrix for the trim frame is not crucial, it turns 
out that the choice of MRP to parameterize attitude de- 
viations from trim is crucial in the generation of simple 
control laws with respect to attitude error feedback [6], 
as will be evident later. The kinematical equations gov- 
erning MRP are 

e = - [(1 — e T e)/ 3 x 3 — 2e x + 2ee T ] u> Fh / Fi (6) 

The variable uj Fb / Ft denotes the angular velocity of Fb 
with respect to F t expressed in Fb frame. Hence if 
the trim frame does not rotate with respect to F a then 
LO Fb / Ft = uj, but in general, u> = lo — ui . Notice 

that the above MRP are directly related to the more fa- 
miliar Euler parameters [5]. Unlike MRP however, the 
Euler Parameters are not independent but their rota- 
tional kinematics are simpler. In summary, the attitude 
error can be conveniently described as || e|| 2 , about a given 
trim attitude. 


Figure 2 shows a flight dynamics and control schematic 
for a general NFTP vehicle with the use of both bias 
momentum wheel and a set of variable speed CMG’s. 
The vector thrusting (and differential throttling) com- 
mands include four fan speeds (w™ d ) and four vane an- 
gles ( Q cmd ) to control platform attitude (/?) , and velocities 
{v and ui) . These vector thrusting commands for the fans 
and vane angles are primarily used for generating trim 
conditions and for open loop command tracking. For the 
general configuration with CMGs, a secondary but never- 
theless crucial function of the vector thrusting commands 
will be to stabilize the CMGs so that their saturation 
tendencies will be mitigated, as described in more detail 
later. To improve performance, fan speed regulators and 
vane servos are used. 

Notice that the Ducted Fan/ Vane aerodynamics and 
gravity directly influence system momenta p 7 h , while the 
CMG Dynamics have no influence on system momenta; 
rather, the CMGs internally redistribute angular momen- 
tum to control the platform angular velocity u> and its at- 
titude (3. A key advantage is that the CMG can generate 
internal control torques to effectively and very reliably 
redistribute angular momentum to control the platform, 
independently of Ducted Fan/Vane aerodynamics. This 
of course mitigates the uncertain but significant effects 
of unsteady aerodynamics induced by rapid control sur- 
face motions for attitude stabilization. This property of 
this control approach will particularly benefit smaller ve- 
hicles during hovering or operations at low airspeeds un- 
der significant wind turbulence because they will respond 
rapidly due to their smaller rotational inertias and will 
therefore require control effectors that can reliably gener- 
ate control forces and moments at higher bandwidth for 
attitude stabilization. 

While the bias momentum wheel is regulated to a cer- 
tain constant speed, to provide directional stability 
in open loop, the CMG subsystem will consist of closed 
loop torque commands to gimbal and wheel motors for 
variable speed CMGs. These commands will be generated 
by a control law which will integrate the CMGs with the 
vector thrusting control system for a given momentum 
wheel augmented platform. 


3 Controller Synthesis 

An important consideration in the selection of a control 
law design approach in this study is global asymptotic 
stability (GAS) due to the inherent complexity in the 
nonlinear, multibody, dynamical system. In addition, 
the large excursions in the gimbal angles, significant gy- 
roscopic coupling, and large attitude angle motions lim- 
its the reliability of system stability, analysis, and design 
methods for linear control systems. This is of particu- 
lar concern when the vehicle needs to cope with strong 
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Figure 2: Flight dynamics and control schematic. 


turbulent winds while hovering or flight at low airspeeds, 
typical of VTOL flight operations. For these reasons, 
we consider using Lyapunov’s direct method in the con- 
trol law synthesis to satisfy GAS. However, in practice 
even guaranteeing GAS is only a necessary condition for 
a successsful flight since controller effectiveness and sat- 
uration can still lead to loss of vehicle. Of course in lieu 
of an expensive flight test program, this dilemma can be 
addressed, with some confidence, only through extensive 
simulations based on specific detailed flight hardware, for 
specific flight environments. 

In this paper, we formulate our control problem along 
the lines of [6], [7], [5], and [8]. In particular, we take ad- 
vantage of the control law simplifications through the use 
of MRP to parameterize the attitude excursions about 
trim, and using the particular form of the excursion error 
as suggested in [6]. For the CMG portion of the con- 
troller, wheel and gimbal torques are commanded, rather 
than the CMG output torques. In proceeding this way, 
we follow the approach in [9], whereby the CMG dynam- 
ics appear explicitly in the control law design. In doing 
so, we avoid separate vehicle control and CMG steering 
laws; instead, the joint control law continuously keeps the 
CMGs desaturated while stabilizing the airframe. This 
continuous desaturation is important, because a vehicle 


of this type operates in an environment of pervasive and 
unpredictable disturbances from atmospheric, manuever, 
and other causes, and cannot safely anticipate opportu- 
nities for dedicated “momentum dumping” events of the 
type typically seen when using CMGs in spacecraft con- 
trol. 

Define the normalized state variables 

Y = D~ 4 X (7) 

where Y £ 1Z 14x1 , and a normalizing non-singular diago- 
nal matrix is given by 

D = diag(ti n ,o; n ,? 7 „,fl„) (8) 

Consider the following Lyapunov function 

V A \y t WY + 2k o ln(l + e T e) (9) 

where k 0 > 0 is a scalar weight specifying the importance 
of vehicle attitude excursions relative to the remaining 
state variables, and {W = W T £ lZ 14xl4 ,W > 0} is 
a symmetric positive definite matrix which defines the 
significance of individual normalized states in the control 
problem. Using the motion equations in 2 and vehicle 
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attitude kinematics in 6, after some algebra, we can show 
that 


V = Y t { W(MD)-\Z + U) + D 


/ 03x1 

k Q e 
04xi 
\ O4X1 


( 10 ) 


One simple approach to enforce Lyapunov stability con- 
dition on the closed loop system is to choose the control 
input U £ 7Z 14x1 such that 


V = -Y t PY, P = P T > 0 (11) 


which requires the following nonlinear state feedback con- 
trol 


U = -MDW~ 1 [PY + D 


( 03xi 

k Q e 

04xi 

\ 04 X i 


\ 


- z 


( 12 ) 


The weight matrix P £ 1Z 14 x 14 can be used as a matrix 
of control gain knobs to optimize during the controller 
design iterations. 

Finally, note from condition 11 that only Lyapunov 
stability is guaranteed since the important state e need 
not be zero for V to vanish. So we state the following 
result that guarantees GAS. 


Theorem 1 The solution X = 0, e = 0, of the closed 
loop system given in equations 2, 6, and 12 is GAS. 


The significance of the above zero solution is that it cor- 
responds to our trim solution. 


Proof of Theorem 1 We use an extension of Lyapunov 
stability theorem found in page 509 of [2]. It says that the 
system is asymptotically stable if there exists a V > 0 such 
that V < 0 and such that there is no non-trivial solution 
such that V = 0 . 

From the choice of V in 9, V > 0, and from equation 
11, V < 0. It is also clear that V = 0 if and only ifY = 0. 
So, what remains is to prove that the set of all solutions 
which also satisfies V = 0 consists of the null solution. 

From equation of motion 2, any solution requiring 
Y = 0 means that 

Z + U = 0 (13) 

since D is a constant, nonsingular matrix. It also follows 
from the feedback equation 12 that 


U = —MDW~ 1 D 


( ° \ 

k Q e 

0 

V 0 


- Z 


(14) 


so that if we substitute 14 into 13, the solution variable e 
must satisfy 


( ° 

MDW~ l D fc ° e 

V 0 


\ 

/ 


— Le — 0i4xi 


(15) 


where L £ 7 ?. 14x3 is a full rank matrix since MDW~ l D £ 
7 ^ 14 x 14 a i SQ a f u n ran h matrix. The only solution then 
is e = 0. The asymptotic stability applies globally since 
the above proof applies independent of initial conditions. 


4 Simulation Results 

The attitude control approach described earlier is applied 
to the vehicle pictured in Figure 3. More details of the 
test vehicle are described in the earlier companion paper 
[1]. The simulations presented below show that globally 



Figure 3: The NASA Flying Test Platform. 

stabilizing nonlinear control laws can be developed which 
are capable of integrating thrust vectoring, bias momen- 
tum wheel, and active use of variable speed CMGs. In 
particular, the controlled performance possible beyond 
bias momentum stabilized system is improved by the use 
of a set of variable speed CMGs, analogous to the case of 
advanced spacecraft control systems in use today. 

In our earlier paper [1] , we demonstrated the fact that 
a system with significant level of bias momentum to pro- 
vide directional stability is significantly more robust than 
a corresponding system with no bias momentum. So, in 
all simulation cases which follow, a significant level of bias 
momentum is assumed. Two sets of control designs are 
compared. The first set consist of controllers that use 
only vanes and fan speeds for active control. The second 
set uses variable speed CMGs in addition to vanes and 
fan speeds for active control. 

4.1 Crosswind and turbulence model 

To simulate wind disturbances during hovering or low air- 
speed operations, crosswinds are assumed which impinge 
on a vertical cylinder attached to the upper part of the 
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platform. The resulting aerodynamic drag on this cylin- 
drical column is used to simulate the crosswind and tur- 
bulence on the vehicle. The resulting disturbance model 
consists of forces and torques in the pitch and roll axis. 
The “strong” winds are intended to simulate “Calm to 
Strong Breeze/Near Gale” and its simulated sample mean 
wind speed is 15 Knots with a standard deviation of 18, 
as seen in Figure 4. 
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Figure 4: Simulated strong wind. 


4.2 Controller Law Synthesis 

Based on the control law described earlier, controllers are 
designed to stabilize and hold the platform attitude un- 
der strong turbulent crosswinds as described previously. 
For a fixed set of performance weights, W, k a , a set of 100 
randomly chosen symmetric, positive definite gain matri- 
ces, Pi 4 X i 4 and Pqx6 are used in a closed loop simulation. 
Only the root mean square (RMS) values corresponding 
to the top 10 cases in terms of requiring the least amount 
of RMS torque from vector thrusting system are shown 
in the following figures. Hence, the control performance 
discussed below is only a small sample of what is possible, 
let alone optimal control. 

Figure 5 shows the tradeoff between the RMS atti- 
tude error and the RMS torque command to the vector 
thrusting/differential thrusting system. This figure shows 
that the performance of systems using CMGs can vary 
widely (if carefully designed at random!) but there exists 
a significant set of controllers using CMGs that requires 
significantly less torque command to the vector thrusting 
system. Figures 6 and 7 demonstrate the same result ex- 
cept in terms of RMS thrust command to vector thrusting 
and RMS angular velocities. 
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Figure 5: Attitude Response vs. Torque Command for 
thrust vectoring. 

4.3 Minimum torque controller 

Figures 8 and 9 show vehicle response comparisons for two 
cases. In both cases, they represent the minumum RMS 
torque (command to the vector thrusting system) confi- 
grations, first case (solid line) respresents the no CMG 
configuration, and second case (dashed line) with CMGs. 
These two cases corresponds to the best control gains in 
figure 5. Figures 8 compares the vehicle response under 
strong turbulent winds. The vehicle attitude excursions 
about trim (shown in the bottom right subplot as e) is 
in terms of MRPs and shows robustness in the attitude 
response. Nice exponential decay is shown for both cases 
when the strong winds disappear at the end of 10 sec- 
onds. The excursions in the gimbal angle rates and CMG 
wheel speeds are very reasonable. Figure 9 compares the 
external net thrust and torque commands to the thrust 
vectoring system (shown as the two subplots on the top 
figure). These are command signals beyond trim for sta- 
bilization. Notice that the CMG case (dashed line) re- 
quires significantly less torque commands to the thrust 
vectoring system than the system without CMGs (solid 
line). The figure also shows a reasonable level of CMG 
wheel and gimbal motor torque commands. 


5 Conclusions 

This paper has examined the use of internal angular 
momentum transfer devices - bias momentum wheel 
and variable-speed CMGs, to enhance stabilization of 
a thrust-levitated VTOL vehicle. A nonlinear control 
law was formulated that satisfies conditions guarantee- 
ing global asymptotic stability, and gains we designed via 
random selection. Nonlinear simulation showed signifi- 
cant enhancement of closed-loop attitude performance in 
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Figure 6: Angular Velocity Response vs. Torque Com- 
mand for thrust vectoring. Figure 7: Attitude Response vs. Thrust Command for 

thrust vectoring. 


hovering when subjected to cross-winds disturbances. It 
was also seen that the CMGs could reduce the amount of 
thrust vectoring actually needed for stabilization. Future 
work will concentrate on techniques for efficient controller 
design, and on formal criteria for sizing the CMGs and 
bias momentum against the thrust vectoring capability 
for this class of vehicle. 
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Figure 8: Vehicle response to strong wind, CMG (dash), No CMG (solid). 
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Figure 9: Control response to strong wind, CMG (dash), No CMG (solid). 
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